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I.  Introduction. 


Much  work  has  been  done  recently  in  the  area  of  the  superposition 
of  independent  stationary  point  processes,  extending  earlier  results 
summarized  in  Cox  and  Lewis  (1966,  Ch.  8)  and  Cinlar  (1972).  Various 
approaches  have  been  utilized  in  deriving  the  probabilistic  properties 
of  the  superposition  process,  including  the  development  of  expressions 
for  the  description  of  counts  of  events  occurring  in  some  interval  and 
for  ;.he  description  of  the  properties  of  intervals  between  events 
starting  at  an  "arbitrary"  event.  Following  Lawrence  (1971)  we  call 
these  the  synchronous  intervals  of  the  process. 

The  present  paper  considers  the  problem  of  describing  the  second- 
order  joint  moment  structure  of  the  synchronous  intervals  for  the  special 
case  of  the  superposition  of  identical  Erlang  renewal  processes. 

Notationally ,  let  p  be  the  number  of  component  processes  to 
be  superposed  and  let  each  component  process  be  an  Erlang  renewal  pro¬ 
cess,  i.e.  have  independently  distributed  intervals  between  events,  each 
interval  having  the  9ame  Erlang  probability  density  function 

,  k  k-1  -kx/p 

g(x)  -  (-)  - — p— -  (x^O).  (1.1) 

The  superposition  of  p  of  these  •'ndependent  processes  gives  the  super¬ 
posed  process  to  be  studied.  We  use  the  description  "Erlang  process" 
rather  than  "Gamma  process"  to  indicate  that  the  shape  paremeter  k 
in  (1.1)  is  au  integer. 

Downton  (1972)  considered  the  problem  of  deriving  the  distribu¬ 
tions  of  synchronous  intervals  in  this  superposition  process  and  utilized 
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both  an  Imbedded  Markov  property  and  An  Integral  equation  approach  to 
develop  limiting  forms  for  the  distributions,  i.e.  for  large  times  of 
observation  and  large  numbers  of  superpositions.  Dovnton  was  concerned 
with  applications  to  reliability  studies;  note  that  the  superposed 
process  Is  the  overall  failure  pattern  for  a  system  consisting  of  p 
renewable  components  operating  in  series,  i.e.  the  system  fails  whenever 
a  component  fails  and  is  renewed. 

Lawrance  (1973),  amongst  other  things,  extended  Downtons  results 
and  showed  that  the  synchronous  intervals  between  events  in  the  super¬ 
posed  process  are  correlated.  He  calculated  the  first  two  serial 
correlation  coefficients  for  specific  cases  of  superposed  Erlang  renewal 
processes,  and  obtained  joint  and  univariate  distributions  of  synchronous 
intervals . 

Barnett  (1970)  used  simulation  to  investigate  some  of  the 
correlational  properties  of  superposed  processes  in  the  sre:ific  area 
of  "studying  the  supply  pattern  of  goods  de^ ivered  to  a  supply  depot 
by  a  fleet  of  vehicles."  His  processes  may  be  considered  to  be  Erlang 
processes  with  large  values  of  k. 

Here  we  develop  somewhat  more  general  expressions  than  those 
gi^en  previously  for  the  serial  correlations  and  the  spectrum  of  syn¬ 
chronous  Intervals  it  superposed  Erlang  renewal  processes.  Specific 
forms  for  the  spectra  are  given  for  a  few  cases.  A  general  expression 
for  the  generating  function  of  the  distribution  of  counts  is  used  to  compute 
the  spectrum  of  intervals  and  rerial  correlations  in  other  cases.  Some  general 
conclusions  about  the  shape  o'  the  spectrum  of  the  intervals  in  the  super¬ 
posed  process  are  given;  these  utilize  results  of  Enns  (1970)  to  the 
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effect  that  If  intervals  in  the  component  processes  have  increasing 
failure  rates  (decreasing  failure  rates) ,  then  so  do  the  marginal  syn¬ 
chronous  intervals  in  the  superposed  process.  Enns  (1970)  aleo  discussed 
serial  correlations  and  tried  to  relate  it  to  the  properties  of  the 
failure  rates  of  the  intervals  in  the  component  processes. 
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2.  Preliminaries . 

Cox  and  Lewis  (1966,  Ch.  4)  used  the  fundamental  Identity  of 
point  processes. 


(N(t)  <  k)  «*  (Sk  >  t)  (k  -  1,2,...)  (2.1) 

and  Palm-Khinchine  formulae  to  relate  the  generating  function  of  the 
asynchronous  counts,  ^(g,t)  “  to  the  spectrum  of  synchronous 

intervals  {X^}  in  a  stationary  point  process.  The  development  leans 
heavily  on  Laplace  transforms;  for  (1.1)  we  have 


g(s) 


g(x)e  Sx  dx 


0- 


(1 t  8  +  *>'*. 


(Re(s)  i  0) 


(2.2) 


which  is  a  rational  function  in  s. 

Note  that  the  parametrization  for  (2.2)  gives  an  expected  value 
of  y;  the  intervals  X  can  be  construed  to  be  the  sum  of  k  independ¬ 
ent  exponent- ially  distributed  variables,  each  with  mean  y/k  -  3. 

Denoting  the  Laplace  transform  of  the  generating  function  for 
counts  in  an  interval  (0,t]  following  an  arbitrary  event  (synchronous 
case)  in  the  point  process  by  4f(€;s),  we  have  for  a  renewal  process 
(Cox,  1962,  p.  37) 


*/ 


JLl  &ll8_ 

s[l-{g*(s)}] 


(2.3) 


This  expression  is  simple  looking  but  in  general  difficult  to  invert. 
Using  9*(£,s)  the  Palm-Khinchine  formulae  (Cox  and  Lewis,  1966,  Ch.  4) 
give  the  Laplace  transform  $  (£;s)  of  the  generating  function  for  the 
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number  of  events  in  (x,x+t),  where  x  is  an  arbitrary  tine  in  a 
stationary  process  (i.e.  the  asynchronous  case).  The  relationship  is 
(Cox  and  Lewis,  1966,  Ch.  4,  p.  68) 

MSjs)  -  ±  +  4*U;s).  (2.4) 


If  p  renewal  processes  are  superposed,  we  need  to  develop  a  general 
expression  for  the  transfora  of  the  generating  function  of  the  asyn- 

*(n) 

chronous  counts  in  the  superposed  process,  denoted  by  +  (£;s). 

This  is  because  the  spectrum  can  be  obtained  as  (Cox  and  Lewis,  1966, 
p.  77) 


f+(w)  -  ~  (1  +  2  l  Pk  cob  (loo)} 
k»l 

V  24* (0,0+)  -  v  ’ 


(2.5) 


where  the  p^'s  are  the  serial  correlation  coefficients  for  synchronous 
intervals  k  lags  apart,  i.e.  corrtX^Xj^} ,  and  the  denominator  is 
var(Xi)/E(X1)  -  var(X1)/y. 

In  principle  it  is  easy,  in  the  case  of  superposition  of  p 
Identical  processes,  to  obtain  +  ^(£;s)  or  f  ^(£;t),  since  in 
the  time  domain 


♦  (p)(c;t)  -  E(C 


N(P>(t)}  _ 


N,(t)+...+N  (t) 
EtC  1  P  } 


P  N  (t) 

-  n  Etc  > 
i-i 


-  t*(e;t))p. 


(2.6) 
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However,  If  $  (£>s)  is  known,  as  for  the  renewal  process,  to  obtain 
*(p) 

♦  r  (£;s)  it  is  necessary  to  perform  the  operations  (using  2.6) 

♦*(e;s)  L->  *(£;t)  -M(p)a;t) 

where  L  1  denotes  inverse  Laplace  transformation  and  L  a  Laplace 

transform.  Obviously  the  procedure  would  be  cumbersome  for  large  p, 

* 

even  if  $  (£;s)  could  be  inverted  to  give  4>(£»t). 

Section  4  of  this  paper  is  concerned  with  the  development  of 

♦  ^(S;t)  or  its  transform  in  the  special  case  of  the  superposition 
of  independent,  identical  Erlang  renewal  processes,  while  Section  5 
develops  the  spectral  density,  f+(w)  for  the  superposed  process. 

Section  6  gives  specific  functional  forms  for  f+  (uj)  and  the 
p^'s  in  the  cases  of  k  ■  2  and  p  -  2,3,4;  the  results  are  surpris¬ 
ingly  simple.  In  Section  7  computational  results  for  f+(w)  are  given, 
using  the  general  results  of  Section  5,  for  k  -  2,3  and  p  *  2, 3,..., 9. 
These  illustrate  the  general  form  of  the  spectrum  of  the  superposed  pro¬ 
cess.  Some  results  for  e><jperpositions  of  non-identical  Erlang  processes 
are  also  given,  corresponding  to  cases  discussed  by  Lawrence  (1973) . 

In  Section  3  we  review  some  results  of  Enns  (1970)  and  show  how 
they  can  be  used  to  derive  some  properties  of  the  spectrum  of  superposed 
renewal  processes,  not  necessarily  Erlang.  All  the  results  are  used  in 
Section  8  in  the  discussion  of  the  general  shape  of  the  spectrum  of 
intervals  of  superposed  Erlang  renewal  processes. 
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3 .  Ena's  results  and  the  initial  point  on  the  spectrum. 

It  is  clear  that  the  marginal  time  between  events  In  a  stationary 
super  pored  process  X^  will  have  mean 

E(X(p))  -  E(X)/p, 

where  E(X)  is  the  common  oea‘.>  time  between  events  in  the  p  component 
processes.  Let  the  spectral  density  (2.5)  of  the  synchronous  intervals 
in  the  superposition  process  be  denoted  by  f£p^  (to)  and  the  coefficient 
of  variation  squared  by  C2U(P))  “  var(X(p))/{E(X(p)}2.  Then  we  have 
the  following  relationship  between  the  asymptotic  slope  of  the  variance 
time  curve  V^(t)  •  var{N^P\t)}  and  the  initial  point  on  the 
spectrum  (Cox  and  Lewis,  1966,  p.  78): 


lim 

t-K» 


dt 


IT 


C2(x(p))  (p) 

E(X(P))  + 


(W-). 


(3.1) 


Now  V^(t)  *  pV(t),  where  V(t)  is  the  variance  time  curve  for  the 
component  processes,  and  if  the  component  processes  are  renewal  pro¬ 
cesses  it  is  well  known  (see  Cox  and  Lewis,  1966,  p.  31)  that 


lim  V' (t)  -  lim 
t-x»  t-*®» 


dYlil 

dt 


C2(X) 
E(X)  ‘ 


Putting  these  results  together  we  get,  for  superposed  renewal  processes, 


P 


V'(«)E(X^) 
ir  C  (X^) 


.  P.  Cjtt)  x  ._EjX)  .  -Si SSL- 

*  E(X)  p  C2(X(P))  Tt  C2(X(P)) 


(3.2) 


There  are  two  consequences  of  this  result. 
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Enns  (1970)  showed  that  If  X  has  an  Increasing  hazard  rate 
distribution  (IHR) ,  then  so  does  ,  and  if  X  has  a  decreasing 
failure  rate  distribution  (DFR) ,  so  does  X^ .  Thus,  using  results 
of  Barlow  and  Proschan  (1964)  we  have  that 


C*(x*p))  il  if  X 

is 

IHR, 

C2(X(p))  *  1  if  X 

is 

DFR. 

Consequently 

„  C2(X<P>)  +  * 

if 

X  is  IHR, 

if 

X  is  DHR. 

.  C2(X<P))  + 

Note  that  for  the  Erlang  density  function 

(1.1) 

,  X  is  IHR  if  k  >  1 

and  DFR  if  k  <  1,  and  C2(X)  -  When 

k  - 

1  w?.  have  an  exponential 

density  and  a  Poisson  process  for  the  component  processes  and  the  super¬ 
position  process;  the  spectrum  then  has  the  constant  value  1/ir  and 
f+(Of)  -  f|p)«*)  -  C2(X)/w  -  1/ir. 

The  second  consequence  of  (3.2)  is  that  in  the  limit  a&  p  -*■  «, 
when  the  superposed  process  goes  to  a  Poisson  process  and  C  (X^)  1, 

the  spectrum  of  Intervals  is  tied  down  at  w  ■  0: 

limit  f£p)  (0+)  -  . 

p-H» 

For  the  Erlang  process  this  value  is  1/ (ik) .  Thus  the  spectrum,  which 
is  converging  to  the  flat  spectrum  of  the  limiting  Poisson  process,  has 
a  "notch"  in  it  at  w  near  zero.  This  notch  corresponds  to  the  fact 
that  if  we  could  observe  a  superposed  process  well  past  the  mean  of 
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intervale  in  the  component  processes,  the  variance  tins  carve  would  pot 
be  equal  to  t/E(X) ,  i.e.  would  depart  fro*  linearity.  In  other  worda 
low  frequency  effecto  aud  long  tine  effecta  correspond. 

It  ia  not  known  if  f^(Ot)  goes  ncnotonically  fro*  1/v  at 
p  *  1  to  l/(«k)  as  p  ■*■  «.  Computational  result  a  are  given  in 
Section  7.  It  would  also  be  interesting  to  know  the  derivative  of 
f£p^(tt)  at  at  ■  Of. 


v.T 
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4.  General  results  for  the  r.aneriting  function  for  an  Erlang  renewal 
process . 

For  the  Erlang  renewal  process,  using  (2.2)  and  (2.3)  and  letting 
X  *  k/y  -  1/0,  we  have 


♦f <C;«) 


l  -  (-^-)k 

2Xfsi 

■a-e<xs» 


(X+s)k  -  lk 
a{(X+s)k  -  Xk}  * 


(4.1) 


Hence,  using  (2.4), 


♦*<?;•) 


ks 


(“H(X+s)k 


,  (X+s)k  -  Xk__ 
s{(X+s)k  -  CXk> 

-  ?Xk}  +  {-^~^H(X+s)k  -  xk} 

(X+s)k  -  {Xk 


(4.2) 


(4.3) 


where  P(s)  and  Q(s)  are  polynomials  in  s. 

By  expanding  and  collecting  terma  in  P(s)  one  finds  that 

P(s)  -  -  Q(s)  +  (Q(»)  +  Xk(C-l)>  (4.4) 

-  aQ  +  a1s  +  a2s2  +  ...  +  a^a1*"2  +  a^s**1, 

so  that  P(s)  is  a  proper  polynomial  in  s  of  order  k  -  1,  with 


aQ  -  Xk_1{k  +  (C-l)(^)), 
ax  -  Xk”2{^}{k  +  (C-l)^)>, 
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a2  -  X 


k-3,<k-V 


+  (C-i)(^;’, 


*k-2  “  A{k  +  ^T}  » 


The  first  derivative  of  Q(s)  with  respect  to  s  Is 


Q'(s)  -  k(X+s)k_1. 


(A.  5) 


Nov  let  {t. . Tj,}  be  the  k  distinct,  complex-valued  roots 


€  *  i«e« , 


Tj  ■  C  (J  -  l,...,k). 


The  roots  s^  of  Q(s),  i.e.  *  0,  are  clearly 


8j  -  X(Tj-l)  (J  -  1,2 . k). 


(A.  6) 


and  these  too,  like  ,  are  distinct  and  complex  valued. 

We  are  nov  In  a  position  to  use  a  standard  Inversion  formula 
(Gardner  and  Barnes,  19A2,  p.  338),  for  the  ratio  of  two  polynomials 
when  the  denominator  has  distinct  roots: 

1  k  P(s.)  s.t 


L  lQ(»r  ^0*0^7 


We  note  that  from  (A. A)  and  (A. 6) 

.k-1 

p/8  )  .  h —  (£iL.\2 
riV  it 

and  from  (A.S) 


(J  “  1 , • • • »k)  , 


Q'CSj)  -  MX+Sj) 


k-1  kX 


(j  -  l,...,k). 


(A.  7) 
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Consequently,  the  probability  generating  function  for  the  asynchronous 
number  of  counts,  N(t),  in  an  Erlang  renewal  process  Is 

♦(C;t)  -  L“1{^*(5;») ) 


k  P(s.)  s  t 

m  T  -  ■pJL»  ■  a  J 

-W-^T^TF*  1  •  (4'8> 


We  note  too  that  the  expansion  for  the  variance-tin*  curve  of  an 
Erlang  renewal  process  obtained  by  Serfllng  (3070)  is  obtainable  from 
(4.8). 
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5 .  General  Result*  for  the  Spectrum  of  Counts  for  Pooled  Independent 
Erlang  Renewal  Proceaaea. 

We  consider  here  the  caae  In  which  each  of  the  p  component 
Erlang  renewal  process  have  the  sane  scale  parameter  X  and  shape 
parameter  k.  The  probability  generating  function  for  counts  in  the 
pooled  process  starting  at  an  arbitrary  time  is,  from  (2.6), 

♦  (p'a,t)  -  (♦ce,t)}p. 


Let  J  -  •  •  »Jp}»  be  an  index  set.  Using  the  result  (4.8)  from 

Section  4, 


♦  <l,)  (5»t)  - 


TT^F  •“V“| 


W  Mi  ^1 


X(  l  (T  -I)i 

.  jei  J 


(5.1) 


Taking  the  Laplace  transform  of  $^(C,t)  we  have 


k 

•  I 


Ll 


j  -1  {  8  -  X  l  (T  -1) 

p  V  j€J  3 


(5.2) 


Vie  note  the  following  in  (7  >5): 


(ii)  the  denominator  the  second  expression  of  (2.5)  for  f+(u>) 
is  constant  with  respect  to  «. 
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Consequently,  1£  we  ean  evaluate  ^p\e**,Of)  +  4*^p\e-lw,0+)  for 
0  *  u  s  *,  and  Integrate  it,  then  (i)  will  give  us  the  value  of 
*{24*Cp\o,0+)  -  E(X) }.  Thus  ve  will  have  f|p)(o>)  for  0  *  «  *  ir, 
and  also  var(X'*').  Mow  when  a  -  Of, 


**(p)  a.o*-) 


.ijigaif  {  ...  sj,L2^!j 


(5.3) 


We  are  Interested  in  the  cases  g  -  e*“  and  g  »  e”*10.  We  note  that 


-  c  -  2  +  r1. 


(5.4) 


Thus, 


u=m 


-I£zii2 

lid  g 


g-e 


lb) 


-lb) 


g-e 


.^-e-2  +  e-  -2  (1-cos  u) 


Recall  that  is  a  k^*1  root  of  g.  We  will  now  use  the 


notation 


,  th  ,  id 

Tj  •  a  k  root  of  e  , 

t ^ 1  *  a  kC^  root  of  e  *w. 

In  particular  we  assign  indices  to  these  roots  in  the  following  manner; 


J  *  lg.ee ,kg 

(5.5) 

j  *  1  |  e  e  e  )ke 

(5.6) 

Thu.,  Tj'  -  Tj-1. 


15 


Taking  real  and  Imaginary  parts, 


Re(tj)  -  Re(tj  1 )  -  , 

K^)  -  -Ktj')  -  8in(^J&)  . 


(5.8) 


(5.9) 


Because  Tj '  "  Tj  »  w  aee  th^t 


■  ~;~2  +  v1  • 


Thus ,  for  any  Index  set  J, 


n  (_  >i)2  "  n  (x  *  _i) 2  "  H  _i 

J6I  UJ  *  j€J  Uj  '  j6J  ^  -  2  +t  *} 


j6J  -2{  1-cob  (2rJ+a>)  ) 

We  now  consider  the  terms  l/][(l-x^)  and  1/^(1-t^') 

In  (5.3)  when  £  is  replaced  by  e ^  and  e  respectively. 


I  u-V 

J€J  ] 


l  {1-Re(r  )}  -  1  l  I(x  ) 

J€J  3  J&  3 

][  {1-Re(x  )}+  1  l  I(t.) 

_i£ _ i _ m _ 1 _ 

|  |  l  {1-Re(x .) }  -  1  l  I(x  )||2 
J6J  3  J6J  3 


(5.11) 


I  (1-Ti)  l  { 1-Re(x, ')}  -  i  I  I(r') 
j€J  J  j€J  3  J€J  3 

l  (1-Re(x,)}  +1  [  I(x  ') 

_  1€J _ 3  1€J  3 

||  l  { 1-Re(x, ’)}  -  1  l  I(x.')||2 
J€J  3  j€J  3 
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l  {1-Re(0}-1  l  I(tJ 

JM  J  Jei  3 

|  l  {l-Re(T.)}-i  J  1(0  || 
J6J  3  j€J  3 


(5.12) 


Thua ,  for  *ay  index  set  J, 


I  (1-T)  I  (1-t  *) 
j€J  2  j€J  2 


2a1 


{5.13} 


where 


and 


aT  -  l  {l-Re(tj) } 


j€J 


-  I  (1-cos  (^)>, 
J6J 

b.  -  l  Kx  ) 

3  J€J  2 

.  I  . 

J€J 


Let 


c  -  n  ( 1-cos  (2,rj+u)) )  . 
J  J6J  C 


(5.14) 


Gathering  results  (5.3),  (5.4),  (5.10),  (5.13)  and  (5.14)  we  have 


f<p)(»>  ~  ♦*(p)(.1‘>)  +  ♦*(P)(.-1“,W) 


■  <!-«=»•  “>P  I  •••  I  (aW)c  '  (5- 

Alt  P  Jx-1  j  -1  '8FDJ3CJ 


15) 


Because  Sj,  l .  and  Cj  depend  only  on  the  indices  in  J  and 
not  on  their  order,  the  number  of  distinct  terms  to  be  computed  for  each 
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value  of  «  is  equal  to  ,  the  figurate  number  for  k  types 

of  objects  taken  In  sets  of  alee  p.  The  conputatlon  tine  for  a^, 
bj  and  Cj  la  proportional  to  p,  so  total  conputatlon  time  should 
be  approximately  proportional  to  p  times  the  figurate  number. 

A  FORTRAN  program  was  written  to  compute  fj^  (w)  at  Intervals 
of  G.05ir.  Various  values  of  k  and  p  having  figurate  numbers  no 
greater  than  56  were  used.  Computing  times  were  less  than  2  seconds 
per  run  on  the  IBM  360/67  computer.  Some  results  of  the  computations 
will  be  shown  in  Section  7.  First  we  obtain  some  specific  results 
for  small  p  and  k. 
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6.  Specific* raaulta  for  k»  2  and  p  ■  2,3,4. 

The  generating  function  ♦  ?^(C;t)  for  the  superposition  of  p 
identical  Erlang  renewal  processes  is  analytically  tractable  for  k  -  2 
and  small  p. 

From  (4.8),  with  a  ■  -X(l+5^2)  and  b  ■  -X(l-5^2>»  we 

obtain 

♦  tt;t)  -  — Y72  {-(1“51/2)2  eat  +  (1+51/2)2  ebt).  (6.1) 

45  ' 

By  taking  {♦(5*t)}*>  we  obtain  the  generating  functions  for  the 
superpositions  of  p  such  Erlang  renewal  processes. 


(i)  For  the  case  p  -  2,  we  have 

(♦(C.t)}2  -  ♦<2)«,t)  -  ^  {(1-51/2)4  e2at 

+  (1+£1/2)2  ebt  -  2 (1-5) 2  e(a+b)t},  (6.2) 


whose  Laplace  transform  is 


U;a) 


_i_  I  U-e1/2)l,|  an1'2)*  .  2(i-e>2} 

165  (  s  -  2a  s  -  2b  s-a-b  j  * 


(6.3) 


From  (6.3)  we  can  obtain  the  spectrum  of  intervals  using  (2.5).  The 
result  is 


<2>(»)  -  ±  (i  -  ^p) 


(0  £  to  £  r) . 


(  v.4) 


From  (2.5)  we  have  immediately  that  ■  -1/10,  *  0  for  k  >  1, 

and  the  spectrum  is  that  of  a  first  order  moving  average  process.  This 
generalises  results  of  Lawrance  (1973),  who  gave  p^  -  -1/10,  but  no 
higher  correlations. 
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(11)  For  p  ■  3  the  generating  function  of  the  counting  function  of 
the  superposition  process  is 

♦  <3)(C;t)  - 575—  {-(l-e1/2)6  e3at  +  3(l-C1/2)2(l-5)2  e(2a+l,*t 

<4C1/Z) 

-  0(l+e1/2)2(l-O2  e(**2b)t  +  (i+51/2)6  e3bt>,  (6.5) 


and  its  transform  is 


♦*(3)«;s)  - 


1  J-U-E 


3/2)  s  -  3a 


l~S1/2)2(1-P2 

a  -  2a  -  b 


s  -  *  -  2b 


(l+£1/2)6j 

•  3*  J* 


(6.6) 


Again  using  (2.5)  we  obtain  for  the  spectrum  of  intervals 


~(3)/  v  9  9123  -  48  cos  w  -  9  cos  x  ,, 

f+  <“>  ‘  25i\ - U  -  9  cob  hj - }  »“*’)•  <6'7> 

This  is  the  spectrum  of  a  mixed  second-order  moving  average  and 

first-order  autoregressive  process  for  which  the  serial  correlations 

are  expressed  by  three  constants,  since  for  k  i  2  p,  ..  -  p8,  . 

Rr  l  K 

k 

Equivalently,  for  k  i  2,  -  C8  . 

In  the  present  case  p^  -  -0.1044447,  p2  -  -.0316033, 

P3  -  -0.0035118,  or  C  -  -2.5596b  and  8-1/9. 


os2  ft)  I 

“  f 


(ill)  Skipping  the  details,  we  get  for  p  -  4  the  result  for  the 
spectrum 

*(4)/  ,v  8  9871  -  502  cos  u  -  57  cos2  id  -  24  cos3  u>  (  ^  _  _x  /4. 

f+  <“>  ■  379VI - i~'--a~co."» - }  <°  *  “  *  *)  (6-8> 
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This  Is  the  spectrum  of  mixed  third -order  movlrg  average  and  first- 
order  autoregressive  process,  so  that  the  serial  correlations  may  be 
expressed  by  four  constants,  rince  for  k  >  3,  ■  gp^.  Equivalently, 

for  k  st  3,  p^  -  Cgk. 

For  the  case  p  -  A,  k  -  2  we  get 

Px  -  -0.0979540  ; 

P2  -  -0.0442821  ; 

P3  -  -0.0150263  ;  C  -  -.961683  ; 

P4  -  -0.0037597  ;  8  -  1/4  . 


In  Figure  1  we  show  the  spectra  f^p)  (u>)  for  the  superposition 
of  Independent  renewal  processes  with  Erlang  (k~2)  interval  distribu¬ 
tions.  The  calculations  for  p  -  2,3,4,  where  p  is  the  number  of 
processes  superimposed,  are  exact  results  using  (6.4),  (6.7)  and  (6.9) 
respectively.  The  result  for  p  ■  5  was  obtained  from  the  general 
results  of  Section  5. 

Note  that  the  initial  points  of  the  spectra  seem  to  be  decreasing 
monotonically  with  p  to  the  value  'few)  *  generally  foT;  the 

other  a)  the  spectrum  seems  to  be  converging  toward  the  flat  spectrum 
of  the  Poisson  limit,  i.e.  — . 
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7.  Discussion  of  Results . 

The  results  (6.4),  (6.7)  end  (6,8)  giving  a  simple  structure  to 
the  spectra  of  the  superposed  processes  for  p  -  2  and  k  »  2,3  and 
4  are  at  first  quite  surprising.  However,  some  such  structure  could 
be  anticipated  from  (S.1S)  which  gives  the  spectrum  f^(u)  as  a 
rational  function  In  cos  u.  It  has  not  besn  possible,  as  yet,  to 
relate  the  parameters  p  and  k  to  the  order  of  the  moving  average 
and  autoregressive  structure,  or  to  give  simple  expressions  for  the 
spectrum  or  serial,  correlations.  Knowing  the  order  of  the  structure 
of  the  spectrum  would  be  a  help  in  knowing  how  many  serial  correlations 
one  woii.d  have  '  •>  compute  from  (5.15)  for  specific  p  and  k.  The 
correlations  are  computed  by  calculating  f+(u>)  at  a  sufficiently  fine 
grid  and  then  inverting  the  series  using  a  FFT  algorithm.  (See  Cooley, 
L.  wis,  Welch,  1970.) 

There  are  other  theoretical  questions  raised  by  the  results  which 
will  be  discussed  else’ here.  Briefly  though,  Downton  (1971)  discussed 
the  possibility  of  looking  at  this  particular  superposition  process  as 
a  Markovian  progression  through  pk  parallel  channels,  each  channel 
having  k  stages.  This  may  explain  the  structure.  For  p  ■  k  *  2  it 
is  In  fact  possible  to  show  that  the  Intervals  are  not  only  uncorrelated 
but  Independent  for  lags  greater  tiian  2. 

We  discuss  now  the  computational  results. 

In  Figure  1  the  spectrum  is  given  for  the  superposition  of  Erlang 
(k"2)  processes  for  several  values  of  p,  the  number  processes  super¬ 
posed.  This  was  discussed  In  the  previous  section.  A  similar  graph 
Is  given  for  k  «  3  in  Figure  2.  Again  as  p  increases  the  initial 
point  of  the  spectrum  decrease  is  converging  slowly  to  the  value  1/v 
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but  the  remainder  of  the  spectrum  which  would  obtain  for  all  u  for 
the  limiting  Poisson  process. 

Figure  3  by  contrast  has,  for  fixed  number  p  ■  2  of  superposed 
processes,  graphs  for  different  values  k.  As  k  increases  and  the 
component  processes  become  more  regular  with  C2(X)  decreasing,  there 
is  a  larger  peak  at  ir,  or  equivalently  at  a  period  of  T  ■  2.  Thus 
the  alternation  of  intervals  induced  by  the  superposition  becomes 
pronounced  in  the  spectrum  of  Intervals. 

In  Figure  4  the  case  p  »  3  is  examined  for  different  values 
of  k.  The  tripling  due  to  the  superposition  produces  a  more  and  more 
marked  peak  at  u  ■  2ir/3  (period  T  ■  3)  as  the  component  processes 
become  more  regular. 

Formulae  such  as  (5.15)  can  be  developed  for  non-identical  component 
processes.  Tn  Figures  5,  6  and  7  several  cases  of  the  pooling  of  a 
Poisson  and  Erlanglan  process  are  given;  the  marginal  olstributlons 
were  given  by  Lavrance  (1973).  These  are  important  as  models  for  the 
case  in  which  there  is  one  dominant,  fairly  regular,  process  corrupting 
the  superposition  of  many  sparse  component  processes  which  sum  (almost) 
to  a  Poisson  process.  Note  carefully  the  vertical  scale;  departure 
from  a  flat  spectrum  is  small  and  would  be  hard  to  detect  iu  real  data 
unless  the  Gamma  process  were  very  regular.  The  spectrum  of  counts  would 
probably  be  much  more  informative  in  this  case. 

A  short  table  of  computed  coefficients  of  variation  of  the 
marginal  intervals  X  in  the  superposed  process  is  given  in  Table  1. 

These  came  out  of  the  computations  of  the  spectrum. 
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POOLED  EBLAHG  PROCESSES 

L  '  1  ' 

Coefficient  of  Variation 

r  ‘  “  -  "i  i  'I  ■»  f  i  ‘ 

p  -  2  j  ;  3  4  5  1  6  7  8  9 

to*2 

3 

4 

5 

6 

7 

8 
9 


.791 , 

.833 

.860 

*79 

.893 

;904 

.913 

.920 

.714 

.779 

.819 

.847 

.867 

.882 

.894 

-V." 

.675 

.754 

El 

fJJ 

WM 

) 

Hi 

.652 

<740 

.  t , 

•  ■ 

i  r 

HI 

.637 

.732 

f  r.  •  J,**' 

V 

'  i'-i'. t* 

.627 

.620 

.615 

TABLE  1 


8.  Conclusion. 


*-*Vr'r.  ■'  •■■*'• ' *'  ■ 

The  fact  (Cox  and  Lewis,  1966,  Ch.  8)  that  the  variance-time 

t  t  '  1  X  i  V  .1'’  .  J  ^  #<* 

curve  and  spectrum  of  counts  in  the  superposition  of  Identical  processes 
is  a  scaled  version  of  the  corresponding  components  is  very  useful  if 

...  -V-  ! 

P  is  known.  The  crux  jof  moat  analysis  though  is  that  p  is  not  known 

,  ;  .  .  ’  ••• . f  ■  ■  •  '  .  ; 

and  must  be  determined.  The  spectra  of  intervals  computed  here  should 

.  I  ’  4-  t  ~  —  ' 

be  a  very  useful  tool  in  helping  to  determine  p,  and  in  verifying 

s  1  ■  >  i 

i  .  .*■  < 

assumptions  such  as  those  made  in  Cox  and  Lewis  (1966,  Ch.  8)  in 
analyzing  nerve-pulse  data. 
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CAPTIONS 

Figure  I.  Computed  spectra  of  intervals  toe  superposed  Erlang  renewal 
processes  (k«2)  for  several  values  of  p,  where  p  is  the  number 

of  independent  processes  superposed.  For  p  ■  2,  3,  and  4  these  are 
exact  results  using  (6.4),  (6.7)  and  (6.9). 

Figure  2.  Computed  spectra  of  intervals  for  superposed  Erlang  renewal 
processes  (k»3)  for  several  values  of  p,  where  p  is  the  number 

of  independent  processes  superposed.  For  all  p  the  results  are 
obtained  from  the  computational  formulae  of  Section  5. 

Figure  3.  Here  the  number  p  of  superposed  processes  is  held  fixed  at 

2,  while  the  parameter  k  in  the  pooled  ITrlang  processes  is  increased. 

As  k  increases  the  Intervals  in  the  individual  processes  become  more 

regular  and  the  peak  in  the  spectrum  of  the  superposition  process  at 

u  "  n  becomes  more  pronounced. 

Figure  4.  Here  the  number  p  of  superposed  processes  is  held  fixed  at 

3,  while  the  parameter  k  in  the  pooled  Erlang  processes  is  increased. 

As  k  increases  the  intervals  in  the  Individual  processes  become  more 

regular  and  the  peak  in  the  spectrum  of  the  superposition  process  at 

u  -  2w/3  becomes  more  pronounced. 

Figure  5.  The  spectrum  of  intervals  of  superposed  Poisson  and  Erlang 
renewal  processes  is  plotted  for  Increasing  values  of  the  parameter  k 
in  the  Erlang  process.  The  Poisson  and  Erlang  process  mean  values  are 
held  constant. 

Figure  6.  The  spectrum  of  Intervals  of  superposed  Poisson  and  Erlang 
renewal  processes  is  plotted  for  increasing  values  of  the  parameter  k 
in  the  Erlang  process.  The  Poisson  and  Erlang  process  mean  v  ilues  are 
held  constant. 

Figure  7.  The  spectrum  of  Intervals  of  superposed  Poisson  and  Erlang 
renewal  processes  is  plotted  for  increasing  values  of  the  parameter  k 
in  the  Erlang  process.  The  Poisson  and  Erlang  procecs  mean  values  are 
held  constant. 
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